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%4 ' We examine analytically the ghost propagator Dyson-Schwinger Equation (DSE) in 

the deep IR regime and prove that a finite ghost dressing function at vanishing momentum 
is an alternative solution (solution II) to the usually assumed divergent one (solution I). 
We furthermore find that the Slavnov- Taylor identities discriminate between these two 
classes of solutions and strongly support the solution II. The latter turns out to be also 
preferred by lattice simulations within numerical uncertainties. 
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How to deal with the IR behaviour of QCD ? There are three main types of approach: 



\^ | • Dyson Schwinger equations (DSE) and especially the untruncated one concerning the 

ghost propagator. 

(N 

^ . • Ward-Slavnov- Taylor identities (WSTI) 

o , 

00 , • Lattice QCD simulations (LQCD). 

o 

Until a few years ago, there was a clear contradiction between the standard DSE solution 
and LQCD results. If we call F(q 2 ) (G(q 2 )) the ghost (gluon) dressing function, the standard 
DSE solution (later labelled as solution I) predicts that F 2 (q 2 )G(q 2 ) goes to a non-vanishing 
constant when q 2 — > (see for instance [1] and references therein). LQCD indicates on the 
contrary in an unambiguous way that F 2 {q 2 )G(q 2 ) -> when q 2 -> [2, 3]. The standard 
solution implies [4] also that G(q 2 )/q 2 does not diverge when q 2 — ► while F(q 2 ) diverges 
at least as fast as (q 2 )' 1 ^ 2 - Regarding lattice QCD results, they have long been compatible 
with an IR-diverging F(q 2 ), although definitely at a much slower pace. This discrepancy has 
been tentatively charged to different types of lattice artifacts. However more recent LQCD 
data obtained in large volume simulations [15,16] show that under those conditions the ghost 
dressing function IR exponent of (assuming F(q 2 ) 2 — o (q 2 ) aF ) lies in the vicinity of 0. 

Now, it was proven in [3,5,6] that: 

• there exists a second class of solutions to the DSE (later labelled as solution II) which 
implies that F(q 2 ) goes to a non-vanishing constant when q 2 — ► and does not constrain 
F 2 (q 2 )G(q 2 ). 
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• the WSTI implies under very plausible assumptions that F(q 2 ) goes to a non- vanishing 
constant when q 2 — > 0, which imposes the solution II of DSE. 

Thus, the convergence of the three methods towards a finite non-vanishing ghost 
dressing function is very impressive. 

Furthermore, a recent numerical study of the DSE using the LQCD gluon input finds that 
both cases of solutions (I and II) are found depending on the strong coupling constant which 
is a free parameter in this exercise [7]. Solutions exist when the coupling constant is smaller 
than (or equal to) a critical value. In the general case the solutions which come out belong 
to type II, but for the critical coupling constant one finds the solution I. It was also proved 
that for an appropriate coupling constant the resulting ghost dressing function (belonging to 
class II) fits very well with lattice results. 

Concerning the gluon propagator the analytic methods are not so constraining. WSTI, 
under a regularity hypothesis for the longitudinal-longitudinal-transverse gluon vertex func- 
tion, predicts a divergent gluon propag ator when q 2 -> [3, 8] while LQCD seems to point 
towards a finite non- vanishing gluon propagator at q 2 = (see, for instance, [10]). A very slow 
divergence of the gluon propagator, not easy to see in LQCD, might solve this discrepancy. 

In this paper we wish to present an analytic study of the ghost propagators of both 
solutions I and II in the deep infrared in the context of the DSE. We also will carefully 
scrutinize the relationship between DSE, WSTI and LQCD solutions. In section 2, the ghost 
propagator DSE is properly renormalised and analysed in the deep IR regime. The two types 
of solutions are obtained in section 3 and their implications put clearly on the table. In 
section 4, we discuss what WSTI tells us and section 5 is devoted to briefly review the LQCD 
results for the ghost propagator. We conclude in section 6. In appendix [X] we show how 
the ghost propagator DSE that we exploit in the next section can be generally inferred from 
WSTI. 



2 The ghost propagator Dyson- Schwinger equation 

We will examine the Dyson-Schwinger equation for the ghost propagator (GPDSE) which 
can be written diagrammatically as 



q-k 




i.e., denoting by (resp. G^) the full ghost (resp. gluon) propagator, 



(^ (2) )^) = S ab k 2 (1) 



- glials j ^Ff(g)(^)f^(-9,fe;9-^)(^)(G( 2 ))g(g-^ 

where T stands for the bare ghost-gluon vertex, 



Ff c (-q, k-q-k) = ig f abc qlr u , u (-q, k;q-k) 

= ig f abc ( q^q, k) + (q - k) u H 2 {q, k) ) , (2) 

where q and k are respectively the outgoing and incoming ghost momenta and go is the 
bare coupling constant. Let us now consider eq. (pQ) at small momenta k. After applying 
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the decomposition for the ghost-gluon vertex in eq. ([2]), omitting colour indices and dividing 
both sides by k 2 , it reads 



F(k 2 ) 



1 + 9qN c 



d\ I F(q 2 )G((q - kf 



(2nY 



\q-kf 



(k ■ qf 
k 2 



Hi(q,k) 



(3) 



It should be noticed that, because of the transversality condition, H2 defined in eq. ([2]) does 
not contribute for the GPDSE in the Landau gauge. 

2.1 Renormalization of the Dyson-Schwinger equation 

The integral equation eq. ([3]) is written in terms of bare Green functions. It is actually 
meaningless unless one specifies some appropriate UV-cutoff A and performs the replacements 
F(k 2 ) — > F(k 2 ,A) . . . It can be cast into a renormalized form by dealing properly with UV 
divergencies, i.e. 



G R {k 2 ^ 2 ) 
F R (k 2 ,fi 2 ) 



Z-V,A) 5o 2 (A) 
Z^(fj, 2 ,A)G(k 2 ,A) 
Z^(fi 2 ,A)F(k 2 ,A) , 



(4) 



where fi 2 is the renormalization momentum and Z g , Z3 and Z3 the renormalization con- 
stants for the coupling constant, the gluon and the ghost respectively. Z g is related to 
the ghost-gluon vertex renormalization constant (defined by Tr = Z\Tb) through Z g = 
Zi(Z^ 2 Z^)^ 1 . Then Taylor's well-known non-renormalization theorem, which states that 
Hi(q,0) + H2(q,0) = 1 in Landau gauge and to any perturbative order, can be invoked to 
conclude that Z\ is finite. We recall that the renormalization point is arbitrary, except for the 
special value fj, = which cannot be chosen without a loss of generality (see, in this respect, 
the discussion in ref [9]). Thus, 



F R (k 2 ^ 2 ) 



Z 3 (^,A) + N c Zx g 2 R {^) E i? (^,/i 2 ;A) 



(5) 



where 



Z R (k 2 ,^ 2 ;A) = J 



<? 2 < A2 d A c 



F R (q 2 ,^)G R ((q-k) 2 ,^) 
q 2 {q-kf 



(fc • qf 

k 2 



Hi, R {q,k;ix 2 



(6) 



One should notice that the UV cut-off, A, is still required as an upper integration bound in 
eq. © since the integral is UV-divergent, behaving as f dq 2 / q 2 {l+\las / (2tt) log (<?/m))) -35//44 - 
In fact, the cut-off dependence this induces in Yi R cancels El against the one of Z3 in the r.h.s. 
of eq. ©, in accordance with the fact that the l.h.s. does not depend on A. 

Now, we will apply a MOM renormalization prescription. This means that all the Green 
functions take their tree- level value at the renormalization point and thus: 



F R {»\f) 



Gr^ 2 ^ 2 ) = 1 



(7) 



2 One can easily check that Z 3 1 (fi 2 , A)E_R(fc 2 , /1 2 ; A) approaches some finite limit as A — > 00 since the ghost 
and gluon propagator anomalous dimensions and the beta function verify the relation 27 + 7 + (5 — 0. 



3 



In the following, H\(q,k) will be approximated by a constant!! with respect to both mo- 
menta and, provided that Hi(q,0) = 1 at tree-level, our MOM prescription implies that 
Hi t n(k, q; /j, 2 ) = 1 and Z\ is a constant in terms of \i. 

2.2 A subtracted Dyson-Schwinger equation 

The renormalized GPDSE, eq. should be carefully analysed. We aim to study the 
infrared behaviour of its solutions and therefore focus our analysis on the momentum region, 
k -C Aqcd> where the IR behaviour of the dressing functions (presumably in powers of the 
momentum) is supposed to hold. One cannot forget, though, that the UV cut-off dependences 
in both sides of eq. © match only in virtue of the previously mentionned relation between 
the ghost and gluon propagator anomalous dimension and the beta function. 

However, in order not to have to deal with the UV cut-off, we prefer to approach the 
study of the GPDSE in the following manner: we consider eq. ([!]) for two different scales, \k 
and Xnk (with n < 1 some fixed number and A an extra parameter that we shall ultimately 
let go to 0) and subtract them 

2,2 ,.2. n vw\2„2;„2 



F R (\ 2 k 2 , /J, 2 ) F R (\ 2 K 2 k 2 ,fi 2 ) 



N C gW ) Zi ^R(X z k z ,fi z ; oo) - Z R (\ z K z k z , oo 



Then the integral in the r.h.s. is UV-safe, thanks to the subtraction, and the limit A — > oo 
can be explicitely taken, 

v a 2,2 2 \ v «2 2, 2 2 v f A / F(q 2 ,n 2 ) ( (fc • q) 2 2 \ 
£ K (A k ,/x ;oo) - Y> R (\ k k ,n ;oo) = J j—^ \ ~ 2 I— p q I 

L (9 -A*) 1 '\ 

An accurate analysis of eq. ([8]) requires, in addition, to cut the integration domain of eq. (J9j) 
into two pieces by introducing some new scale q\ (qo, typically of the order of Aqcd, is a 
momentum scale below which the deep IR power behaviour is a good approximation), 

S it (A 2 fc 2 ,/x 2 ;oo)-S R (A 2 K 2 A : 2 ,^ 2 ;oo) = J m (A) + Juv(A) (10) 

where Iir represents the integral in eq. Q over q 2 < qfi and Itjv over q 2 > q$. Only the 
dependence on A is written explicitly because we shall let it go to zero with k, k and /i 2 kept 
fixed. The relevance of the q^ scale stems from the drastic difference between the IR and UV 
behaviours of the integrand. In particular, for (Xk) 2 <C q$, the following infrared power laws, 

Gm({q - Xk) 2 ,fi z ) = B(^)({q-\kf) aG , (11) 

will be applied for both dressing functions in irR. 
Now, Iir is infrared convergent if : 

olf > —2 IR convergence at q 2 = 

ac > IR convergence at (q — k) 2 = and (q — nk) 2 = (12) 

3 This approximation is very usually used to solve GPDSE. Notice that few lattice data are available for the 
ghost-gluon vertex. However, in a recent computation [2] of that vertex for a zero gluon momentum, Hi (q, q) 
appears to be approximatively constant with respect to q. Of course, more data for different kinematical 
configurations should be welcome to check that approximation. 




We shall suppose in the following that these conditions are verified. We then obtain, per- 
forming the change of variable q — > Xq : 



/ir(A) 



2\ {a-F+a-G) 



((<i-k) 2 y 



(2tt)* 
((,-^) 2 ) QG - 2 



A: 2 



(13) 



The point we have to keep in mind is the fact that the upper bound of the integral goes to 
infinity when A — > 0. This potentially induces a dependence on A whose interplay with the 
behaviour explicitly shown in (I13p we must check. In this limit, the convergence of the integral 
depends on the asymptotic behaviour of the whole integrand for large q. In particular, the 
leading contribution of the square bracket in eq. (|13p behaves as 



[(k - q) 



21 a G-2 



[{nk - q) 



21 a G~ 



~ (q 2 r G ~ 2 (ao-2)(l-«) 



(14) 



A: 



+ + 



A: 2 , Jq- k) 2 



where we expand up to (A; 2 /g 2 )-terms because those in q ■ k, being odd under q^ — ► — g^, 
give a null contribution under the angular integration in eq. (|13f) . Thus, provided that the 
conditions (fT2"]l are satisfied so that Itr, is convergent when q — ► (or q — > A;), its asymptotics 
for small A is 



/go/A 



2(a f +a G )-3 



(15) 



where e, the lower limit of the integral, is a small cut-off that avoids any possible spureous 
singularity appearing after expanding in eq. (fTB"]) for large q. Thus, if ap + a G < 1, the 
asymptotic behaviour of /ir. is given by the power on A in front of the integral in eq. (|15p . 
since the integral itself will remain finite when A — > 0. If op + a G = 1, the integral diverges 
logarithmically as A vanishes. Otherwise, one can change the integration variable back, 
q — > (//A, to get the leading power on A multiplying again a finite integral on the momentum 

q, 



Im(A) ~ < 



' \ 2 ( a G+ a F) if a G + OF < 1 

^ 9 , 



A 2 In A if a G + «f = 1 
A 2 if a G + olf > 1 



(16) 



We have assumed that -f/i is constant when varying all the momenta but (|16j) remains true 
if one only assumes that H\ behaves "regularly" for q 2 ,k 2 < q$ (i.e. is free of singularities 
or, at least, of any singularity worse than logarithmic). 

Let us now consider Jtjv- Its dependence on A, which is explicit in the factor inside the 
square bracket of eq. ([!]), should clearly be even in A : any odd power of A would imply an odd 
power of q ■ k whose angular integral is zero. Since the integrand is identically zero at A = 
and the integral is ultraviolet convergent, it is proportional to A 2 (unless some accidental 
cancellation forces it to behave as an even higher power of A). Thus, in all the cases, the 
leading behaviour of Itr. + Iw, as A vanishes, is given by Iir in eq. (fT6|) . The subtracted 
renormalised GPDSE reads for a G + «f < 1 as: 



F H (A 2 A; 2 ^ 2 ) F R (X 2 ^,^) 



N C g 2 R (p 2 ) Zi /m(A) 



(17) 



for small A. 
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2.3 The ghost-loop integral 

A more quantitative analysis than the one presented in the preceding section can be done 
if we compute exactly the integral /ir(A), defined in eq. (fl"3j) . which gives the contribution of 
the ghost loop to the renormalised GPDSE eq. (fT7|) . If a F + «g < 1 it is possible to perform 
analytically the integral and to find a compact expression for it. In this case, one can write 

l m (\) ~ A(v 2 )B(n 2 ) (x 2 ) {aF+aG) ^(k;a F ,a G )-HKk;a F ,a G ^j (18) 
where A(fi 2 ) and B(fi 2 ) were defined in eq. and 

*(*! *F, ao) = J ^ (g 8 )-'" 1 ((g - fc) 2 p" 2 - g 2 ) , (19) 

provided that op, ag) is not singular to let the subtraction inside the bracket and the 
integral operator in eq. (I18|) commute with each other. Then, following [11], we define 



/M) = TT^L /T^Wfe-*) 2 ) 6 



{k 2 ) 2 + a + b J {2-kY 
r(2 + o)r(2 + fe)r(-q - 6 - 2) 
r(-a)r(-6)r(4 + a + 6) 



(20) 



and obtain 



where 



nJl\a F +ac 

$>(k;a F ,a G ) = — ^2 — H a F,a G ) (21) 



4>(a F ,a G ) = -- (f(a F ,a G - 2) + f(a F ,a G - 1) + f(a F - l,a G - 1)) 

+ -(f(a F -l,a G -2) + f(a F -l,a G ) + f(a F + l,a G -2)) . (22) 
Thus, if a F + a G < 1, 

Iir(A) ~ ^ff^ A 2 )^"^ - k 2 ^+^) 0(a Fl a G ) . (23) 



On the other hand, we know from eq. (|16f) that /ir diverges logarithmically as A goes to 
zero if a F + = 1. In fact, since eq. f|23|) is a reliable result for any a F + < 1 how 
close it may be to 1, such a divergence appears as a pole of a Gamma function of 4>(a F , a G ) 
in eq. (|2ip . We will now compute the leading asymptotic behavior of Ifr as A — > when 
a F + ot G = 1. 

In that case, after performing in eq. (|18p the expansion eq. (|14p and neglecting the term 
odd in g M — » — g M , one finds for the leading contribution 

x J dO sin 4 6 (a G -2 + 2(a G -3) (a G -2) cos y 

~ k 2 (l - n 2 ) A{ ^ B f ] a G (a G - 2)A 2 In A . (24) 

We do not specify the lower bound of the integral over q in eq. (I24p because it necessarily 
contributes as a subleading term, once the ghost-loop integral is required to be IR safe. 
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Finally, if oiF + ac > 1, the leading contribution for 7ir(A) as A vanishes can be computed 
after performing back the change of integration variable, q — > q/X, in eq. (|13p . The first even 
term in eq. (114ft dominates again the expansion after integration, but now it does not diverge. 
Then, if we procceed as we did in eq. (|24p . we obtain 

'»w - *oW> ^ - * 2 > • < 25 > 

for small A and ac + aj? > 1. It should be noticed that Jir in eq. (|25l) depends on the 
additional scale qo introduced in eq. (fTU]) to separate IR and UV integration domains. In 
fact, if one takes qo — > oo, Iir diverges. This means that, when of + «g > 1, the behaviour 
of the IR power laws hampers their use for all momenta in the integral. The finiteness of 
the ghost-loop integral of the subtracted GPDSE can only be recovered after taking into 
account the UV logarithmic behaviour for large-momenta dressing functions @. Furthermore, 
I\jy, also behaving as A 2 , should be also added in r.h.s. of eq. (fT7|) in order to write the 
renormalised GPDSE. Thus, the dependence on A but not the factor in front of it can be 
inferred from the GPDSE with only the information of the assymptotics for small-momentum 
dressing functions. 

3 The infrared analysis of GPDSE solutions 

The starting point for the infrared analysis will be the eq. (I17p for small A, where we will 
try to make the dependences on k, k and A of the two sides match each other. 

3.1 The case ocf ^ (solution I) 

We will first study the case oif ^ 0. Then, the l.h.s. of eq. (fTTl) can be expanded for 
small A as 

(1 (26) 



F R (X 2 k 2 , fj, 2 ) F R (X 2 K 2 k 2 ,fi 2 ) ~ v ; A{pt 

and we will obtain from eq. (|17|) : 



N C g 2 R (v 2 ) Z X A^ 2 ) — ~ 1 , (27) 

where the dependences on k, k and A of the numerator and the denominator should cancel 
against each other. Using for Ijr the form given in after eq. (|16|) . we find three possible 
situations: 



If ac + oiF > 1, applying eq. (|25|) in eq. (|27j) . we are led to the conclusion that only 
ctF = — 1 (and ac > 2) satisfies this last equation and could be an IR solution for 
GPDSE. However, such a solution appears to be in a clearcut contradiction with the 
current lattice simulations. 

If ac + oiF = 1, there is no possible solution because the logarithmic behaviour of Iir 
in eq. (|24l) cannot be compensated by the powerlike one in the denominator of eq. (|2T 



4 The scale qo being of the order of Aqcd, power laws for olf + ac > 1 cannot be then solutions of the 
GPDSE in the MR truncation scheme corresponding to Aqcd — > oo (see, for instance, [11]). The same 



argument holds also for ap + ac = 1, because the ghost- loop integral in eq. (I24|l diverges as A — > for any qo 
fixed as well as for go — * °o for any fixed A. 
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• If otG+otF < 1, eq. (|23j) combined with eq. ([271) implies the familiar relation 2a f + ol g = 
and we have then: 

Nc9W) Z ^iyU {-^a) * 1. (28) 

An immediate consequence of this last condition is the freezing of the running coupling 
constant at small momentum. If the renormalization point, fi, is arbitrarily chosen to be very 
small in order that the dressing functions observe the power laws at k 2 = fj?, one obtains 
A(fi 2 ) = fi- 2a F and B(fi 2 ) = /^ 2q g. Eq. then reads 

N c gW) Z X 4>(-^, aG ) ~ 16vr 2 , (29) 

and should be satisfied for any small value of /U. Consequently, it should remain exact as 
fj, — > and provides the small-momentum limit of the running coupling (which is independent 
of the infrared constants for ghost and gluon dressing functions). 
In particular, if a G = 1, one has 4>(— 1/2, 1) = 8/5 and thus 

N c gUii 2 ) Zi ~ 10^ 2 , (30) 

3.2 The case op = (solution II) 

The case of = is particular in that the leading contributions to the two occurrences 
of F in the l.h.s of eq. (|17f) cancel against each other. We have then to go one step further, 
taking into account the subleading terms. Defining Fir by means of Fm(q 2 , /j, 2 ) = A{fi 2 ) + 
Fm{q 2 ,^ 2 ) we rewrite the l.h.s of eq. © as -(Fi R (X 2 k 2 , fi 2 ) - Fi R (\ 2 K 2 k 2 ,n 2 ))/A 2 (^ 2 ) and 
use the known IR behaviour of /ir(A) from eq !16|) ) in the r.h.s. of eq. (|17|) to get 



A{^ 2 ) + A 2 {^ 2 )q 2 \uq 2 if a G = 1 
A{ji 2 ) + A 2 {n 2 )(q 2 ) a F otherwise 



Fm{q 2 ^ 2 ) = { t/ 2 2 2 (2) „ . (31) 



Furthermore, not only the subleading functional behaviour of the dressing function can be 
constrained but also the coefficient A 2 in eq. (j3Tj) . In fact, if we plug this equation into the 
l.h.s. of eq. (|17|) and expand we obtain : 

(A(u 2 )) 2 , „ ~ ( fc 2 (l - K 2 )A 2 lnA 2 ifo G = l 

A 2 (fj, z ) [ (A z /r) a f (1 — k f ) otherwise, 

Let us consider now in more detail the three possible cases. 

(2) 

• If «g < 1, we obtain from eqs. (I23|32p that a F = a G • Then, 

- N c 9W) Zi m a G ) * Wn 2 , (33) 

where, according to eqs. (|20|22p 0(0, a G ) is given by 

0( °' aG) = 2a G (a G + l)(a G + 2)(l-a G ) {U) 

• Similarly if a G = 1, eq. ([MP applied to eq. (j32l) leads to 



iV G g R (^) Z x ~ 64tt 2 . (35) 
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(2) 



1 



i.e., a ghost dressing function 



• At last, if «g > 1, eqs. (f25j) and ([32]) imply: a 

which behaves quadratically for small momenta, In this case, however, as already said 
the ghost loop cannot be evaluated using the IR power laws over the whole integration 
range and it is therefore not possible to solve the GPDSE consistently, nor even to 
determine the small-momentum behaviour of the dressing functions, without matching 
appropriately those power laws to the UV perturbative formulas. Thus, we are not able 
to derive a constraint for the next-to- leading coefficient, ^(a* 2 )- 

In summary, the GPDSE admits IR solutions with ap = and any a G > 0, provided 



that 



2,1 , . „, .. v'rO' 2 ! A 2, ..2 



A{fi 2 ) ( 1-0(0, a G ) 

5V) 



A ^{ 1+ U, 2 
{ A{^) + A 2 {^)q 2 



16vr 2 

A 2 (fi 2 )B(^ 2 ) q 2 lnq 2 



■ 2 i/A/r)l ( / 2 )"' ; J a G <l 

a G = 1 06) 
a G > I 



where ]j 2 {^ 2 ) = N G g 2 R {^ 2 )Z\ and cj)(0, a G ) is given in eq. The gluon dressing function 

is supposed to behave as indicated in eq. ([TT]). In particular for a G = 1, the gluon prop- 
agator takes a finite (and non-zero) value at zero momentum, B(/jl 2 ), after applying MOM 
renormalisation prescription at q 2 = fi 2 . 



4 The ghost-gluon and three-gluon Ward-Slavnov- Taylor iden- 
tity 

In the previous section, we have analysed the infrared behaviour of GPDSE solutions and 
found that the ghost dressing function can either diverge at vanishing momentum (op = 
—a G /2 with a G > 0) or give a finite value (of = with any a G > 0). As appendix IA1 shows . 
the GPDSE can be derived from the general Ward-Slavnov- Taylor equation [12]. We will 
now invoke the Ward-Slavnov- Taylor identity (WSTI) for general covariant gauges relating 
the 3-gluon, T\ pi/ (p,q,r), and ghost-gluon vertices, 

P\T\^(p,q,r) = ^J-{S pv r 2 -r p r u )f pix (r,p;q) 

(<W<? 2 - q P qfi)T P v(q,p;r) . 



F( P 2 ) 
G(q 2 ) 



to shed some light on that matter [8]. Using for the ghost-gluon vertex the general decom- 
positior@ [13] 

F„p(p,q;r) = S u/X a(p,q;r) - r v q p b(p, q;r) + p u r p c(p,q;r) 

+ r u p p d(p, q; r) + p v p p e(p, q; r) , (38) 



and multiplying by r v both sides of eq. (|3Tj> . one obtains: 

r v pxTx^(p,q,r) = 2 X(q,p;r) [(q ■ r)q p - q 2 r p ] ; (39) 

where 

X(q,p;r) = a(q,p;r) — (r-p)b(q,p;r) + (r-q)d(q,p;r). (40) 



3 We work of course on the energy-momentum shell, so that the relation p + q + r = holds 
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Since the vertex function, T, in the l.h.s. of eq. (I39|) is antisymmetric under p <-> r and A <-> u, 
one can then conclude that [8, 14]: 

F(p 2 )X(q,p;r) = F(r 2 )X(q,r;p) . (41) 

This last result is a compatibility condition required for the WSTI to be satisfied that does not 
involve the 3-gluon vertex and implies a strong correlation between the infrared behaviours of 
the ghost-gluon vertex and the ghost propagator. Now, under the only additional hypothesis 
that those scalars of the ghost-gluon vertex decomposition in eq. (|38p contributing to the 
scalar function X defined in eq. (|40p are regular^] when one of their arguments goes to zero 
while the others are kept non-vanishing, one can consider the small p limit in eq. (|4ip and 
obtain: 

F(p 2 )X(q, 0; -q) = F(q 2 )X(q, -q; 0) + 0(p 2 ) (42) 

This has to be true for any value of q, which implies F(p 2 ) goes to some finite and non- 
zero value when p goes to zero, since neither X(q, 0; —q) nor X(q, —q; 0) are presumably 
zero for all values of q. Rephrased in terms of infrared exponents, the latter argument implies 
that of = 0. 

To reach the above conclusions we did not appeal to the properties of the 3-gluon vertex, 
apart from the symmetry under the exchange of gluon legs. If one assumes in addition that the 
longitudinal part of the 3-gluon vertex also behaves regularly when anyone of its arguments 
goes to 0, the others being kept non-vanishing, a divergent gluon propagator at vanishing 
momentum will be implied [3,6,8]. Of course, as far as it involves a vertex with longitudinal 
gluons which have not been very extensively studied, this last conclusion is not as clean as 
the previous one about the ghost dressing (according to authors of ref. [4] a soft kinematical 
singularity appears for the landau-gauge 3-gluon vertex, however it does not concern our 
proof relying on the regularity of the longitudinal- longitudinal-transverse 3-gluon vertex) . 

In ref. [8], we showed that only a very mild divergence, for example of logarithmic type, 
could be compatible with current LQCD results for the gluon propagator. The IR analysis 
of the previous section can be straightforwardly extended to this case by generalizing 

G m {q 2 ,„ 2 ) = B{f){q 2 ) aG \og"(±>\ , (43) 



the effect of which is to modify eq. f|36|) with 



A{n 2 ) (l - 0(0, a G )°^A 2 (v 2 )B^ 2 )(q 2 re log" (g~ 2 )) a G < 1 

A(» 2 ) (l - { Xl)L 2 A V)B(M 3 ) 9 2 W' +1) (T 2 )) a G = 1 ( 44 ) 
[ A^ 2 )+A 2 ^ 2 )q 2 \og»{q- 2 ) a G >l 

where only the power of the logarithm is then modified. 

Sticking now to the case where olf is zero (for the reasons explained above) and a G is 1 
(as suggested by the lattice results) we are left with 

*fai(«V) = ^R(OV) (1- ^ + ^2 ^m(0^ 2 )W) g 2 log ( - +1) (^)) , (45) 

according to whether there are logarithmic corrections to the gluon propagator [y ^ 0) or 
not {y = 0). Here, M is some scale which is out of the scope of the IR analysis we performed 
in the previous section and, if v = 0, B(/j, 2 ) = G|^(0,/i 2 ) is the gluon propagator at zero 
momentum. 

6 Note also that, for our purposes, it will actually be enough to restrict, and not forbid, the possible presence 
of singularities in the scalar coefficient functions provided that they could be compensated by kinematical 
zeroes stemming from the tensors. 
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5 Ghost propagator from LQCD 



The theoretical study by Zwanziger [17] of the Faddeev-Popov operator on the lattice 
in Landau gauge triggered the first Lattice simulation of the ghost propagator [18] in SU(2) 
and SU(3) gauge theories and the subsequent activity which, mainly for technical reasons, 
was mostly dedicated to the SU(2) lattice gauge theory in the infrared region. It was only in 
the last few years that several studies of the SU(3) ghost propagator focused on its infrared 
region and the Gribov copy problem [2] or on their perturbative [19, 20] and OPE non- 
perturbative [21] descriptions. An unambiguous consensus from LQCD, after all this work, 
pointed that F 2 (q 2 )G(q 2 ) — » when q — > (see, for instance, [2,3]) and, consequently, that 
the solution I is excluded provided that the finite-volume artefacts are indeed under control. 
As a matter of the fact, finite- volume lattice simulations all agree on a ghost propagator 
pretty close to that at tree- level (of — 0) and a gluon propagator not far from being a 
constant at vanishing momentum (ac — 1). 

Very recentely [15,16], simulations on large volumes lattices (with a fair control over the 
finite- volume lattice artifacts) yielded solutions for the ghost dressing function confirming that 
of is indeed in the vicinity of 0. Let us now briefly comment about the ghost propagator 
results from these two papers: 

• The authors of ref. [15] simulated the ghost propagator in 56 4 ,64 4 ,72 , 80 4 volumes 
with an impressive control of the finite-size effects over a huge momentum range from 
q 2 ~ 0.01 GeV 2 to q 2 ~ 10 GeV 2 . They fit an IR exponent, a F = -0.174, that 
appears to be in the vicinity of zero (but negative) and at least much larger that 
the most frequently advocated value (~ —0.5). The fit is however delicate because 
the power behavior is dominant, if ever, only on a very small momentum domain 0. 
Indeed, it is adviseable to try a fitting function inspired from eq. (I44p . Moreover, the 
numerical solution (type II), obtained in ref. [7], after a rescaling because of the MOM 
renormalisation, describes strikingly well the lattice ghost propagator data from ref. [15] 
over a large momentum window, from 0.05 GeV to 3 GeV. 

• The authors of ref. [16] computed an IR ghost propagator exponent, aa{= —of), for 
several 3-dimensional and 4-dimensional lattice volumes (ranging from 140 3 to 320 3 and 
from 48 4 to 128 4 ) and collected the results in their table 1. The values of from that 
table are not only in the vicinity of zero (although being negative) but they approach 
systematically zero when the volume increases. They fit the power behaviour on a small 
domain with two or four momentum data. 

In ref. [7], we showed that the k 2 log(k 2 ) term given by eq. (|45p describes very well the 
behaviour of a numerical solution of the GPDSE, eq. (J3]), for = 1.5 GeV) = 29 (such 
a value corresponds to the best description of our ghost propagator lattice data) and with 
a gluon dressing function taken from a lattice simulation (see Fig. 1 of ref. [7]). We showed 
at the same time that including a logarithmic divergence changes appreciably neither the 
deduced ghost propagator nor the conclusions about the infrared solutions. 

In this same work we analysed in detail the behaviour of the numerical solutions of the 
GPDSE as functions of g 2 (/u = 1.5 GeV) and discovered that a singular solution, behaving 
as 1/q for small momentum (as the relation 2ap + «g = requires), appeared only for the 

specific value g 2 ^ = 1.5 GeV) = 33.198 This solution belongs evidently to what is 

referred to above as class I, with = —1/2 and does satisfy the relation 2cxf + olq = 0. 
Furthermore, the closer g 2 ^ 2 ) to this critical value, the smaller the region near q = where 
eq. (f4"5j) is valid. 

7 The fitted IR exponent is unstable, lying more and more in the vicinity of zero as the momentum domain 
becomes smaller (see Fig. 2 of ref. [15]). 
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6 Discussion and Conclusions 



Thus the present analytical considerations and the previous numerical study con- 
verge towards a consistent description of the set of solutions of the ghost Dyson-Schwinger 
equations: 

• A class of solutions where the ghost dressing function is finite and non zero at q 2 = 
(i.e. cxf = 0), depending continuously on the coupling constant (or equivalently on 
-F(O)). Those solutions do not fulfill the relation 2aF + «g = but appear, for an 
appropriate value of the coupling, to be in very good agreement with the lattice results. 

• An exceptional solution, obtained for a critical value of the coupling is IR-divergent 
with ctF = —1/2. Contrary to the previous ones it satisfies 2a f + «g = but is in clear 
disagreement with the lattice data over a large range of momenta. 

We have demonstrated that the discrepancy between LQCD results (implying unambigu- 
ously that ap — and «g — 1) and the usual DSE solutions (2a f + «g = 0) can be solved 
if the second type (II) of solutions is considered. The existence of this second class besides 
the usual solution (type I) has been proven after carefully renormalising the GPDSE and 
applying a substraction procedure to deal with the remaining (after renormalisation) UV 
singularity. This new solution yields a finite ghost dressing function at vanishing momentum 
while F 2 (q 2 )G(q 2 ) goes to a zero when q — > contrary to what occurs with type I. 

For this (type II) solution, an asymptotic formula of the ghost dressing function is obtained 
that only depends on the IR one for the gluon which is taken as an ansatz in this exercise and 
on the renormalized coupling. The numerical analysis of the GPDSE in ref. [7] proves that 
the type II solution exists for any coupling below a given critical value and that it verifies 
the asymptotic formula. 

The WSTI involving the 3-gluon and the ghost-gluon vertex is particularly useful to 
gain some knowledge about the ghost dressing function: by simply assuming the regularity 
of some of the tensorial components of the ghost-gluon vertex, one can conclude that the 
ghost dressing function is finite and non zero at vanishing momentun. Then, WSTI with the 
mentioned regularity assumption will discard the solution of type I. 

Furthermore, LQCD data point to 2chf + olq ~ 1 (certainly larger than 0). Would one 
wish to reconcile these data with type I solution (2a f + olq = 0), very strong finite- volume 
artifacts would be needed. Such a finite-size effect should strengthen the divergence of a 
ghost propagator behaving at finite-volume like at tree-level and damp to zero the gluon 
propagator. This is very doubtful considering that sizeable discrepancies between lattice and 
solution I appear at momenta of the order of ~ 0.3 GeV. On the contrary, very recent LQCD 
data in large volumes [15,16] show a fair stability as the volume increases and, if any, a trend 
towards solution II (of = 0). This is confirmed by the numerical analysis of ref. [7] which 
proves that both type I and II solutions live at infinite volume for different values of the 
coupling constant. 

It is worth also pointing that some attempts to accomodate lattice data within DS coupled 
equations [23] and within the Gribov-Zwanziger approach [24] led to solutions for gluon and 
ghost propagators that behave pretty much like our solution II does. 

Altogether we strongly believe that the question of the ghost propagator behaviour at 
small momentum is essentially solved. The solution type II of GPDSE avoiding the previous 
discrepancies, the three methods (DSE, WSTI and LQCD) strikingly converge to the same 
result: a finite ghost dressing function at vanishing momentum. The case of the 
gluon propagator needs further study. 
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A The Dyson-Schwinger equation as a Ward-Slavnov- Taylor 
identity 

A very general method to derive Ward-Slavnov- Taylor identities consists in taking advan- 
tage of the transformation properties of 



e G(J) 



J V(A)detMexp i J d 4 x (c - ■^■(d (t A^)(d (t A^) + J£A£ 



(46) 



under gauge transformation (cf. [25]). 

M. is the Faddeev- Popov operator and the notation <, >j indicates that the source term 
J has to be kept, although it will eventually be set to (this is denoted in the following by 
the supression of the J subscript). Taking the derivative of the gauge transformed of eq. ([4*6"]) 
with respect to the gauge parameters leads to the general Slavov- Taylor equation 

i < (d,A a p (x)) >j=< j d A yj;(y) D cb (y)F^ ba (y, x) >j . (47) 

F^ ba (y, x) is the ghost propagator and its presence here is simply due to its very definition 
as the inverse of the Faddeev-Popov operator. If one derives eq. (flT]) with respect to J d (z) 
one gets : 



±<(d^(x))A d p (z)>j = <D db (z)F^ ba (z,x)>j 

+ < Jd 4 yJ c ,(y)D cb (y)F^ ba (y,x)A d p (z)>j (48) 

A first consequence of this relation is the triviality of the longitudinal gluon propagator. 
To see this, it suffices to derive both its sides with respect to z p andto set J to zero. The 
result is 



\ < (d,A a p {x))(d p A d (z)) > = <d p D db (z)F^ ba (z,x)> 

= 5 ad 5 A (z-x) (49) 



To derive the second line we have invoked the fact that dpD db (z), the Faddeev-Popov operator, 
is the inverse of the ghost propagator F^ 2 \ Thus, in momentum space, the general form of 
the gluon propagator for an arbitrary covariant gauge reads 



G$ ab (q)=5 ab 



(50) 



Turning back to eq. (fi8l) and setting J go to zero we obtain 

X - < (d,A«(x))A d p (z) >=< D db (z)F^ ba (z, x) > (51) 
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which is nothing else than the GPDSE. Actually its l.h.s. involves only the longitudinal part 
of the gluon propagator, that we have just seeen to be trivial : 



1 



a 



< (d,A a Ax))A a Jz) >=d p n- 1 (x,z) 



(52) 



As for the r.h.s it can be rewritten as : 

< Df (z)F^ ba (z, x) >=< d p F^ da (z, x)>+i< gf deb A e p (z)F^ ba (z, x) > (53) 

The 3-point gluon-ghost Green's function can be expressed in terms of vertex functions 
and propagators through 



G 



(3)fgh 



(p, q, r) 



: / d A xdHd A ze i v x e irz e iqt < A f p (t)F^ h ( 



z,x) > 



(54) 



F(p 2 ) F(r 2 ) 



G{q 2 ) I £ _ '/,,'/- \ A ,, '/,.'/- 



+ a-. 



ff9hf v (p,r;q)(2n) A S A (p + q + r) 



, r y > P v j ...,7-7 

We Fourier transform eq. (f5"Tj) . use equations (i52|) . (f53|) and (eq. (|54|) ) and obtain 



hp 

k 2 



kp 
G{q 2 



F(k 2 )-gf deb r ba J 

q P q 



eba f d A q F(k 2 )F((k + . 



2 \ P v 



(2vr) 4 k 2 (k + q) 2 

q P qy 



+ a 



(q 



2\2 



T v (k, -k — q;i 



(55) 



The usual form is recovered by multiplying with k p and dividing by F(k 2 ), which leads 



to 



F 



i-gf deb f 

G(q 2 ) 



eba f d A q F{{k + qf 



k v 



(2vr) 4 (k + q) 2 

(qk)q u \ . {qk)q v 
+ a- 



{q 2 ) 2 



T v (k, —k — q;i 



(56) 



This is a general result, valid in any covariant gauge. Of course the a depending (longi- 
tudinal) term disappear in Landau gauge. 

T v (k, —k — q; q) is related to the quantity previously introduced in section H] through the 
relation 



T v (k, -k - q; 



-igkfiT^k, -k - q;q) 



and is usually decomposed into T u (k, —k — q;q) = g [k u H\(k, q) + q u H2(k, q)}. 
Inserting this in eq. (!56p and restricting to the Landau gauge case gives 



F- L (k 2 ) = l + g 2 N c 



d\ F({k + qf 
(2tt) 4 (k + q) 2 



G(q 2 ) ( {qk) 2 _ k2 

q 2 V q 2 



H x (k,q) 



(57) 
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